Electron-electron and spin-orbit interactions in armchair graphene ribbons 
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The effects of intrinsic spin-orbit and Coulomb interactions on low-energy properties of finite 
width graphene armchair ribbons are studied by means of a Dirac Hamiltonian. It is shown that 
metallic states subsist in the presence of intrinsic spin-orbit interactions as spin-filtered edge states, 
in contrast with the insulating behavior predicted for graphene planes. A charge-gap opens due to 
Coulomb interactions in neutral ribbons, that vanishes as A ~ with a gapless spin sector. 

Weak intrinsic spin-orbit interactions do not change the insulating behavior. Explicit expressions 
for the width-dependent gap and various correlation functions are presented. 
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Since the discovery of the anomalous quantum Hall 
effect in graphene [l], 0], the single-layer carbon mate- 
rial has received ever-growing attention due to its un- 
usual physical properties and potential technological ap- 
plications. The unsual properties stem from the peculiar 
band-structure of the material, that at low-energies, can 
be described by a Dirac-type Hamiltonian for massless 
electrons 0] . For this reason, much of the phenomena as- 
sociated with QED can be studied for the first time in a 
condensed matter system in a controlled manner, even at 
room temperatures. At the same time, as graphene holds 
the promise to revolutionize the future of electronics, in- 
creasing efforts arc made to obtain samples with tailored 
properties. In this regard, it is remarkable that in the last 
year, two groups have already succeeded in fabricating 
the first set of graphene ribbons with variable width 0, H| 
for systematic studies. As a consequence, it is essential 
to understand the properties of the material in confined 
geometries as controlled production of graphene nanos- 
tructures is becoming a reality. Moreover, it would also 
be interesting to understand the confinement of Dirac 
fermions in the presence of interactions as it may have 
relevant consequences in other areas of physics. 

The purpose of this paper is to provide precise an- 
swers to two key questions that arise when confinement 
effects are important: a) what are the consequences of 
the enhanced electron-electron interaction on transport 
properties of graphene ribbons? and b) are there metal- 
lic graphene nanoribbons in the presence of intrinsic spin- 
orbit (I-SO) interactions? These two issues have received 
much attention in the cases of two-dimensional graphene 
planes @, 0] and rolled graphene (nanotubes) [1, 
but no detailed study exists at present for nanoribbons 
with defined edges. In this work, we predict that small- 
momentum scattering introduced by electron-electron in- 
teractions opens a width-dependent charge gap in the 
spectrum of half-filled ribbons, that vanishes in the limit 
of wide ribbons. We also present evidence for an in- 
cipient magnetic order in narrow ribbons, in agreement 



with recent experimental measurements [l(J, and show 
that these results are minimally affected by a weak I- 
SO interaction. Unlike previous numerical results 
the approach presented here gives insight into the pro- 
cesses that originate the gap, allowing a systematic ac- 
count of the various scattering processes due to inter- 
actions, and provides expressions for several correlation 
functions. Furthermore, we demonstrate that metallic 
states exist in armchair ribbons of special widths even in 
the presence of I-SO interactions and provide expressions 
for the associated wave functions. This is a fundamen- 
tal difference from reported results on graphene sheets 
where insulating behavior was predicted ■ 

A generic graphene ribbon has a combination of two 
types of edge terminations: armchair and zigzag. Both 
edges arc associated with characteristic transport behav- 
ior in tight-binding calculations: zigzag ribbons are pre- 
dicted to be metallic due to a topological edge state; and 
armchair ribbons are predicted to be metallic or insulat- 
ing depending on the ribbon's width 0, [l2|, 13, 14 1. In 



what follows we focus on finite width armchair ribbons. 

We describe a graphene ribbon using an hexago- 
nal Bravais lattice with a unit cell containing two 
carbon atoms A and B. An atom A is connected 
to its nearest neighbor by 6i = a(0, l/v3) f 62 = 
(a/2)(-l,-l/V3),<5 3 = (a/2)(l,-l/V3); with the lat- 
tice unit vectors given by: dl = (l/2)eyfc(<5j — 6k)- The 
simplest model for electrons in graphene is given in terms 
of a nearest-neighbor hopping tight-binding Hamiltonian 
Hq = Y^<ij> tc\cj+h.c; where i,j label position and spin 
degrees of freedom in real space. The resulting band- 
structure contains six equivalent points in the Brillouin 
zone with vanishing density of states, from which only 
two are independent (the rest are obtained from recip- 
rocal lattice vector translations). We choose these two 
independent points to be located at: K, K' = ±|^(1, 0). 
We work with a spin-dependent basis, the pseudo-spin 
basis fyi> a = (u s A ;ug), where the first (second) compo- 
nent represents the amplitude of the wavefunction at a 
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lattice A(B) site, and s labels the spin. In momentum 
space the Hamiltonian takes the form: 
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where 0=1 + 2 cos(q x a/2)e~ l ^ kya / 2 with a spectrum 
E±(k) = ±t\cf>(k)\. At low energies, the Hamiltonian can 
be reduced even further to its Dirac form, by expanding 
4>(k) around the K, K' points: 
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with v = t\/3/2 (we use h = 1) and k x = q x — K. Con- 
sider a ribbon of length L along the y direction with finite 
width W along x. As shown in the wave function 

that vanishes at the boundaries x = 0,W = W + a/2, 
contains two states at {K ± k x ) and (K' ± k x ) where k x 
takes the values (k x — ^ )W = irn. The corresponding 

energy is given by E ~ ±v^Jk 2 + k 2 . For ribbon widths 

W = (3AI + l)a, a set of linear dispersion (E = ±vk y ) 
states appears with wave functions given by: 
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This solution defines left (\Ev =+ ) and right movers 
(^ r= _) for each spin component. At energies below 
the bulk band-gap (of the order of vf/W ~ OAeV for 
W/a = 10), these linear modes dominate the physics and 
they are the only ones considered in the rest of the paper. 

Effect of Coulomb interactions. It is known that 
Coulomb interactions produce drastic changes in the 
ground state of carbon nanotubes resulting in Luttinger 
liquid physics at low-energies [H, HH . We show next that 
they have similar effects for narrow armchair ribbons. 

The unscreened Coulomb interaction in two- 
dimensional graphene is given by U{x,y) = 



y 2 ) where ao ~ a/2 is the radius 



of carbon p z orbitals and we take the value of n ~ 2.45 
appropriate for a graphene sheet on top of an insu- 
lating substrate [9(. The interacting Hamiltonian is 
H int = (1/2) J dyaly'H where: 

n= J2 Vpviy-y')^ s (yW\y'W^y'W p {y) (4) 

pp' ,ss' 

with effective one dimensional potentials given by 
Vp P >(y~y') = a 2 E„.„' U(Ax nn , , Ay pp ,)\(/)(x n )\ 2 \(f){x n ,)\ 2 ; 
and Ax ni „> = a(n-n' + ±(2y/3d-2y' 73d) mod2 ), Ay PiP / = 
y — y' + pddp^p* and p = + (p = — ) labels the sublat- 
tice A (B). An analysis of the scattering processes shows 
that at low energies, small momentum transfer processes 
(long-range part of Coulomb interaction) dominate for 
reasonable ribbon widths (M > 5), similarly to carbon 



nanotubes [8|, l9( . These forward scattering terms include 
inter- (Ay p _p = y — y' + pd) and intra-sublattices pro- 
cesses (Ay PiP = y - y'). In terms of p p (y) = £ s V'p.sV'p.s, 
(electron density at sublattice p), the Hamiltonian reads: 



H = 



V 

d y[jp(y)p(y) 
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-rr^pp{y)p-p(y)]- ( 5 ) 



f EJp,pW and SV 



where p — pa + Pb, V 

2^2m( V P-p(y™) ~ v pAy™)) measures the difference 
between intra- and inter-sublattice interactions. 

To solve the above Hamiltonian we use bosonization 
17| . As usual, bosonic fields are introduced : ip r ^ s = 
'' r,s expl— ikpry— i-v/7r/2($ c +s<f> cr +r0 c +rs0 (T )l where 
77 is a Klein factor, <& c and are charge and spin bosonic 
fields (0 C and Q a are their duals). In terms of the total 
charge density (p = %=■) and the spin current (J = 



f) y e 



2tt 



^=) the Hamiltonian #($, 6) reads: 
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- „ „ / dy[cos(\^$ c + 4k F y) - cos V8Tr$<r} (6) 



where 



c/ c 



= v/K c , a> K c = l/Jl + and K> 



1/ J 1 — . This Hamiltonian describes two decoupled 
sine-Gordon models for charge and spin sectors. Notice 
that within this approximation, metallic ribbons of width 
W = a reduce to quantum spin chains [l7j]- Analysis of 
Eq. © shows that away from half-filling (kpa = tt/2), 
spin and charge sectors are gapless. At half-filling, K a > 
1 and the spin field $ CT remains gapless, but K c < 1 and 
the charge field becomes gapped. The gap A is estimated 
in the self-consistent harmonic approximation [lH ]: 



A = — ' 



a \ 2ttv c 
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Figure Q] shows the dependence of the band-gap and the 
K c , K a parameters on the ribbon width for a ribbon 
length L = 10 6 \/3a ~ 400/ito (the dependence on the 
ribbon length enters in the definitions of V and SV as 
shown above). The band-gap takes sizable values for 
narrow ribbons (A ~ l30meV for a W ~ 5nm ribbon), 
although it decreases dramatically with ribbon width as 
A « l/W a , with a — ► 1 in the infinite length limit. The 
logarithmic dependence of a on L results from the ap- 
proximation of keeping only long-range scattering terms 
and can be dealt with the formalism presented in a sys- 
tematic way. These results are in qualitative agreement 
with numerical calculations [ll| and experimental mea- 
surements 0, i|. Recent DFT calculations [Il[ on arm- 
chair ribbons have found a width dependent gap that 
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FIG. 1: Left: Charge band-gap of an armchair ribbon of 
length L ~ 400^im as function of its width W. The data 
show behavior A w 1/W a with a = 1.6 (see discussion in 
text). Right: K c and K a as function of the ribbon's width. 



vanishes as 1/W in the limit of wide ribbons. The var- 
ious widths considered in that case included insulating 
and metallic ribbons, thus comparison with our results 
applies only for metallic ribbons (see Ref.[ll|, Fig. 2b, 
N=3p+2). Interestingly, although DFT calculations in- 
clude effects of electron-electron interactions, the opening 
of a gap in Ref. was associated with relaxed atomic 
positions at the edge with carbon atoms passivated by 
H atoms. Measurements on lithographically patterned 
graphenc ribbons with different widths have also revealed 
non-metallic behavior with a width-dependent gap, ex- 
plained as a consequence of quantum confinement. Com- 
parisons of these results with theoretical predictions how- 
ever remain at a qualitative level due to the undetermined 
nature of ribbon's edges in the samples used. 

Unlike numerical methods, the approach used pro- 
vides expressions for various correlation functions 
(T t O(t, y)O^(0, 0)), where O is an operator represent- 
ing relative inter-sublattice charge (Oc = Pa — Pb) or 



x, y, z) densities. These have 



27r$ c )cos(v / 27r$<r)) 
27r$ c ) sin(y27r$ CT )) 



spin (psi =S\-Sh,i = 
the following expressions: 



{O c Oc) = N(sm(y' 4 
(Os*O s *) =N(cos(y'- 
(O s ±O s ±) = N(cos(y' + V^c^^ -)) (8) 

with y' = 2kpy and N = At half-filling, the charge 
field is locked at $ c = and the dominant correlation 
functions arc (O s *O s *) ~ y~ K " and (0 S iO S i) ~ y- 1 /** 
(i = x,y). These results suggest enhanced magnetic 
correlations in the direction parallel to the graphenc 
plane for narrow ribbons. Away from half-filling how- 
ever, the dominant correlations depend on the values 
of K c and K a which are length dependent in the for- 
malism presented. For < K c < 1 the most rel- 
evant operators are Oc,Os^>v with correlation func- 
tions decaying algebraically as (OcOc) ~ y^ K "~ K " and 
{Os*,vOs*-y) ~ y- K "- 1 / K "- respectively. In the regime 
K c < 4p the dominant correlation function is Occ = 
(Pa - Pb)(pa - Pb) that decays as (OccOcc) ~ V~ AKa - 
This regime might be achieved only for narrow ribbons. 



Effect of intrinsic spin-orbit (I-SO) interaction. An- 
other interesting characteristic of graphene is that its 
lattice symmetry supports a particular spin-dependent 
second-neighbor hopping interaction, known as the I-SO 
interaction. This interaction is responsible for the ap- 
pearence of spin-filtered edge states 0, [l9[ which have 
received much attention because of their role in spin 
quantum Hall physics and their relevance for spintronic 
circuitry. As proposed in the I-SO interaction is an 
extension of Haldane's original model [2(| for the quan- 
tum Hall effect with total zero magnetic flux that includes 
both spins. Its Hamiltonian representation is given by: 



Hso 



7 cr z 
-7CT 2 



(9) 



where 7 = 2t'(s\n{q x a) — 2 sm(q x a/2) cos(v / 3fc y a/2)), t' is 
the I-SO coupling, and a z is the spin Pauli matrix. An in- 
tuitive picture for this term is shown in Fig. [2j the orbital 
motion of an electron hopping in sublattice A(B) couples 
to the spin of an electron sitting in the enclosed site of 
sublattice B(A). Since the hopping is spin-dependent, 
the interaction preserves time reversal invariance. After 
expanding the Hamiltonian around the (K, K') points, a 
linear dispersion mode along the y-direction is found: 



fyi- = Nsm e 

3a 



Tsio(x-W'/2) 



1 



(10) 



where N = 



270 



is the normalization factor, 



7o 



sinh 70 W 

3\/3t' , and s = ± labels the spin. This expression re- 
duces to the known result in the limit of semi-infinite 
ribbons [2l| and it is in sharp contrast with predicted 
results on graphene sheets [7[. As Eq. (fTU|) indicates, 
these new modes are localized on the edges and are spin- 
filtered states: theprobability currents along the edges 
are spin-polarized [6(. The right panel in Fig. [2] shows an 
example of such state. The velocity of this linear mode 
is equal to the Fermi velocity of the free system since in 
the presence of the I-SO term, k x acquires an imaginary 
part k x = 47r/(3a) -M70 that preserves the linear disper- 
sion relation. As done previously, it is possible to define 
an effective Hamiltonian along the y-dircction and write 
it in terms of new boson fields (& c , $^) and their duals. 
By including only forward scattering terms, the effect of 
Coulomb interactions on these states is obtained from : 



H = H(&,&) 



dy 



V 1 ' 2 - 26V 1 ' 2 
2^ 



d y *' c By& a (11) 



with the new velocities and Luttinger parameters: 

v' c = v/K' c K' c = + [V 1 - 1 - 25V 1 ' 1 - SV)/2ttv 
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FIG. 2: Left: I-SO term with second-neighbor hopping of an 
electron in sub-lattice A, coupled to the spin of an electron 
on B. Eq.© results from adding clockwise motion and ex- 
changing the role of both sublattices. Here Go = k x a, 6± — 
—k x a/2± k y a /2y/3. Right: probability amplitude for a right 
mover with spin-down (dash-line) and with spin-up (solid- 
line) at the edges of the ribbon as a function of ribbon width. 



Here we have defined V Jr 



q\/3 
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= ^E m V^: p k (y m ) and 
V^{y m )) where 



v pp'(y - v') = a 2 u ( Ax nn>, Ay P j/) cosh«) 
cosh(<,)H) 3+ I- 
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with a l n = 2j x„ - j W + As Eq. (JTTJ) shows, 

Coulomb interactions renormalize the bare couplings and 
introduce a new term d y <&' c d y Q' a that destroys spin- 
charge separation. This term respects time reversal in- 
variance and is present in quantum wires with Rashba 
spin-orbit interactions 22|, 1 231 ] . A rotation by an angle r\ 
diagonalizes the quadratic terms of the Hamiltonian: 



tan 2r\ 



A 



A = v(l 



SV V 1 ' 2 
2nv 7r 



25V 
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(14) 



This rotation modifies the last term in Eq. (|11| pro- 
ducing corresponding changes in correlation functions. 



As it has been proposed [3, |2J|, the strength of the 
I-SO interaction is small compared to the bandwidth 
(i ~ 3eV,t' ~ O.OOlmeU), implying that the rotation 
angle r/ <C 1. As a consequence, the strong coupling 
physics of Eq. Jll]) and the conclusions of the previous 
section remain unchanged in the presence of I-SO. 

To summarize, we studied the low-energy physics of 
armchair nanoribbons in the presence of intrinsic spin- 
orbit and electron-electron interactions within a Dirac 
model. We focused on states with linear dispersion near 
the neutrality (Dirac) points, and showed that they per- 
sist in the presence of I-SO interactions, as spin-filtered 
states localized on the ribbon's edges. These results have 



direct testable experimental consequences and are in con- 
trast with predictions of insulating behavior caused by 
this particular interaction in graphene planes. For half- 
filled systems, we showed that small-momentum transfer 
processes, due to the long-range part of Coulomb interac- 
tions, open a charge-gap in the spectrum, while keeping 
the spin sector gapless. The gap A is strongly dependent 
on the ribbon's width and vanishes in the limit of infi- 
nite ribbons as 1/W in qualitative agreement with DFT 
calculations. Unlike numerical methods [ll|, the formal- 
ism used provides clarification on the origin of the gap 
and gives an explicit expression for its dependence on the 
ribbon width as well as expressions for various correla- 
tions functions. Further analysis of the results suggests 
an enhanced magnetic correlation for narrow ribbons in 
qualitative agreement with interpretations proposed for 
recent experimental measurements [Toj . For graphene 
ribbons away from half-filling, charge and spin sectors re- 
main gapless even in the presence of interactions. These 
results also hold for weak I-SO interactions. 
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ral, S.E. Ulloa and H. Fcrtig. This work was partially 
supported by NSF-NIRT No 0304314. 
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